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Abs tract 


In a perturbed periodic classical motion the angle variable may 
be eliminated either by Kruskal's transformation to "nice variables" 
or — if the system is canonical — by the Poincare-Von Zeipel 
method. For systems that possess a Hamiltonian, the present work 
(l) shows that Kruskal's transformation may be made canonical order 
by order, (2) derives a practical formula for achieving this result 
and (3) shows that the two methods are equivalent and may be matched 
ozder by order. 
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INTRODUCTION 


A basic perturbation problem in celestial mechanics and in guiding center 
motion involves a set of n first-order differential equations which can be 
represented vectorially as 

d£/dt = F(x) (l) 

and which has the following properties: 

(1) F(x) depends on a small parameter 1 and may be expanded 

in it: 

F(x) * F (0) (x) + &F (l) (x) + ... (2) 

(2) In the limit £ -*■ 0 ("unperturbed case") the system may be 
solved and its solution is then periodic in time* 

The problem is then to find an approximate solution valid for small but finite € 
.Mid useful for time intervals of the order of £ periods* 

Kruskal ' ^ devised a method for achieving this* which in many ways resembles 

the method of Bogoliubov and Zubarev and of Krylov and Bogoiiu- 

( 2 ) 

bo* . The calculation in this case proceeds in two steps. First* a transformation 
to "intermediate variables" y(x) is performed* such that in the limit £ — 0 * 
y n is an angle variable linear in time* while the remaining (n-l) components 
of y (which we shall collectively denote by y ) are constant. The equations 
according to which y evolves* derived from eq. (l)* then have the form 

dy/dt = H £ k g (k) (y) (3) 

Jc=0 
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In then* equations, 
periodical functions and 


(k) 

appears in the g only as the angle argument of 
g^ has only one non-zero components, the last one^ 



( 0 , 0 , 


0 , 1 ) 


( 4 ) 


Thv» second step involves a near-identity transformation to "nice variables" z 

2. =» y ♦ 21i ^ 3 r) (5) 

k>t ~ 


such that in the new equations of evolution the transformed angle variable z n 
no longer appears on the right -hand side 

d^dt » 2^, (6) 

k»0 


where as before 


and 


£ - ( £ » *n ) 




(7) 


A general recursion scheme for deriving % order by order has been described 
in an article by Stern^ (henceforth referred to as I ) and ve shall follow 

the notation introduced there, which differs slightly from Kruskal* s (a similar 

(5)v 

scheme for the related Krylov-Bogoliubov expansion has been given by Musen' 

After equations (6) have been derived in this fashion, their first (n-l) components 
constitute an autonomous system for deriving the components of « , which can then 

be independently solved. The problem is thus reduced to one with (n-l) variables. 
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If the system furthermore possesses a Hamiltonian, an additional variable 
may be eliminated by deriving a constant J of the motion which Kruskal terms 
the adiabatic invariant. It is defined as 

1 

p.dq = 


j p*(lq/ i)^) d^ 


( 8 ) 


with the integration performed using an arbitrary canonical set (p, q) , 
over a group of points ("ring”) all of which evolve according to (6) and 
possess the Same z but with values of z n that cover the full range for 
that variable (this property, if once established, is maintained throughout 
the ring's evolution in time). 

If ss itself forms a canonical set, with the momentum conjugate to 
Zq , then one may use this set in (tt), leading to 

1 

J = I Z 1 dz n 58 z l ( 9 ) 

0 

which is a great simplification (if this is not the case, one is forced to 
retrace the transformations x ->• y -► z, until some canonical set, with which 
(8) can be evaluated, is reached). Kruskal did not derive nice canonical 
variables, but he showed them to be possible. Specifically, he proved that in 
any "nice" set, the following Poisson bracket relations are always satisfied 


( z-j] 

L*i • j 1 


in 


( 10 ) 


( 11 ) 


.( 6 ) 


Following Kordheim and Fuea , Kruskal then shove that J and Zq may be 
augmented by (n-2) functions of the z± to form a complete canonical set. 
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In what follows we shall assume that the intermediate variables y ^ form 
a canonical set, in which case a near-identity transformation like (5) can, 
in principle, lead to nice canonical variables. This assumption is not unreason- 
able since, if a Hamiltonian for the system is known, such y^ can in general be 
derived by solving the unperturbed motion via the Hamilton-Jacobi method* 

Ve then 

(1) Show that the freedom allowed by Kruskal's method in the derivation 
of each order of eq. (5) is sufficient to assure the canonical 
character of the "nice variables" s^ . 

(2) Derive a method for obtaining such • 

(3) Show that the result is equivalent to what is obtained by conventional 
perturbation methods based on the Hamilton-Jacobi equation, 

THE POSSIBILITY OF STEP-BY-STEP DERIVATION 


If y is canonical 

y =» (p. q) 

then from (10) 


-[[v 'j] - fo- yj) 


( 12 ) 


If (5) is substituted here, the sere-order part cancels identically and the 
expression remaining inside the curly brackets separates into orders of £• and 
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gives 


(T>/Dz n ) 


Ln w 


m=l 


u) 





= o 
(13) 


Thi3 condition is satisfied for any "nice" set z, . If z, is not merely nice 
but also canonical, then the expressions in the curly brackets of both (12) and 
(l?) are not only independent of z but actually vanish. 

Let us now assume that the expansion ( 5 ) has already been derived and brought 
to canonical form, up to and including order (k-l). Then the first (k-l) orders 
of (13) do, in fact, vanish, leaving (after division by ) 


5i k) . yj - [5 j k) . y t ] + £L^ s) . + - 0 

1 3-1 (14) 

v(k) 

As the next step, one may derive £ and thus extend the calculation of z 
one more order. The equation satisfied by is ( I , eq. 15 ) 


^5 (k) /^y n - h u; (y) 


(k} ( 




(15) 


where depends only on lower orders of and , assumed to 

be known at this stage. If denotes the result of averaging over the 


angle variable y 


n 


< * w > 


V k ^(y) 4y n 


( 16 ) 


one gets (see 1, eq. 19) 
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- J ( >«* * /Hy) 

o 

♦ /^ (k) (y) (IV) 

(k) 

Here ft is an arbitrary additive vector independent of y n , allowed 

by the fact that the derivative of £^) which enters (15) is unaffected by 

(k) 

such an addition. The question now arises whether j>y may be selected so 
as to make the expansion (5) canonical to order k . 

If this occurs, then the 0(l) part of (14) must vanish. It helps here 
to introduce the concept of the conjugate vectors y (I, eq. 4) 

y * ( q . -p ) (18) 


Then 




'"v V ( m ) / 7 — 

J £j /0y i 


(19) 


This has the form of a component cf a curl dyadic in y space, which one 
may call the conjugate curl, with component'’ i'^oted by 


(V * 


( 20 ) 


in an article 




Now it has been shown^by Stem^ (henceforth referred to as II ) that the 
general condition for an expansion (5) to be canonical is 


$ 


(m) 


x (m) (g ) 


^^(m) 

( m — 1, 2, ... ) 


( 21 ) 
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where V isa gradient operator in y space, ~)C ' is an arbitrary scalar 

and f (m) (£ ) are vectors of a certain form, depending on orders of £ lower 

(o) 

than the m-th and on their derivatives. Various choices of fr ' are derived 
in II , all of which satisfy the identity (for canonical £^ 8 ^) 

m-1 

(s) >(m-s) 


& * £ tn) >io - - LW'tr’) 

S=1 

.(a) 


V ' 1 


(22) 


Because the ' in this case are known to be canonical up to and including 

order (k-l), eq. (22) will hold for the lowest (k-l) values of m ; however, 
it will also hold for m = k , for even then the orders of £“ appearing in 
the equation are all lower than the k-th. Substitution in (14) then yields 


CV*z n )[vx[£ (k) - f (k) (£)]j 


♦ 0(s ) 


(23) 


Now let the transformation inverse to (5) be 


y z + 21 (24) 

S— 1 

Given the expansion (5), the may easily be derived (see appendix). 

Alternatively, they may be directly obtained from expanding the relation 
between (3) and (6), in a manner similar to what has been done for the 
expansion (5) in I . This is essentially the method of Krylov and Bogoliubov^ , 

(5) 

as expanded by Musen . Indeed, the elimination of the angle variable by 
Kruskal's method so resembles the Kryl o v-Bogo 1 iubo v approach that the two 
ought perhaps to be regarded as one method (which could be called the Krylov- 
Bogo 1 iubo v-Kruskal method; however, Kruskal's work proceeds past the elimination 
to the derivation of adiabatic invariants). 
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Since the inner part of (23) is a function of y , one must transform 

^>/Dz n = Oyi/^^) ^/^y± 

= Zj{^in + 2 J £' a (^^ 8) /3en)J y t (25) 

or 

'7)/'Dz n =» ^)/"^yn + o(£) (26) 

Equation (23) thus becomes 

W'dyn^V X [^ k) - f (k) (Ojj ♦ 0(£) = 0 (27) 

One may now assume that all variables in (2?) are expressed in terns of y ; 

in that case, each order in £ vanishes separately, including the zeroth. 
Because differential operators commute , this means 

? X • ( 75/7) y n ) [ £ (k) - I (k) (5)]j = 0 (28) 

This integrates to 

0/c5y n ) [ £ (k) - f U) ($ )] = vy (29) 

^(k) _ L ( k \$) = VI. ♦ u (z) (30) 

where T is the indefinite integral of 4* and u is an additive function 
independent of z n , allowed by the integration. 

(k) 

The above condition is satisfied by any £ belonging to a nice set 
of variables which is canonical to order (k-l). For instance, defined 

in U7) will correspond to a certain X and to a certain additive function, 
which may be denoted by u : 

g(k) = f U) (| ) + VI + u(z) (31) 
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If one now chooses in (17) 


r- 


. s 


then 


5 


(k) 


Mk) a 

£ - -a 

f^($) ♦ VI 


(32) 


.00 


and by (2l), £ ' satisfies the condition for canonical variables, making z 

canonical to order k . Thus the requirement can be met* 


PRACTICAL CANCNIZATICM 

(k) 

In order to actually derive the "canonizing" ft * one oust first 
investigate the amount of arbitrariness inherent in that vector. Let 


(k) 


| (k) * yu (k) (3) 


extend the canonical properties to 0(£. i.e. let it satisfy 

7X (k) 


* 


/ K (k> (5) 




(k) 


Then if one replaces / w' 7 with 


(33) 


(34) 


yu'^ - + V J(y) 

with an arbitrary function independent of y , one i^nda 


( 35 ) 
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^(k) + yu,. (!£)(}, „ ♦ V( % W * i ) (36) 

(k) 

which still has the required form (2l). The canonising choice of M ' is 
thus arbitrary within the addition of a conjugate gradient of some scalar <£_ 
which does not contain y n ° 

In order to isolate one operates on (34) with the averaging 

(k) 

operator of (16); since /± * independent of the angle variable, it 

equals its own average, giving 


=■ <(f tk) (5) - ! (l0 _> + (vx w ) 

- <ff (k) <p . v**y * y oc k) > (5v) 

The last terra on the right is a conjugate gradient of some function of y’ , 

~(k) 

i *nA it has already been established that such Vectors , when added to yM , 
do not affect canonization* One may thus drop this part find obtain 



(38) 


To evaluate this, £ ^ must be derived from (17) while f^ may be 
obtained by methods given in II • The most general canonizing additive function 
is then 


* <^f^(|) - + v?(y) (59) 


with arbitrary 
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EQUIVAI£NCB TO CANONICAL PERTURBATION THEORY 


A widely used method for solving perturbed periodic canonical systems is 
due to Poincare and Von Zeipel^'^ 8 ^^ and operates in the following 
manner. First, one expresses the Hamiltonian in terns of the solution 

y « ( p> q ) 

of the unperturbed Hamilton-Jacobi equation, with y n an angle variable and 
jj the conjugate action variable. In the absence of "slowly varying" quantities 
the Hamiltonian then assumes the fora 

H y x ♦ ^L, £ k H (k) (y) (40) 

fcsl 

Next a near-identity transformation to new variables 

* = ( b £ ) 

is sought, produced by the generating function 

<T(P, q) = ZL P i ‘k * ZL q) («) 

k*l 

and having the property that the new Hamiltonian H*(£) is independent of the 
transformed angle variable , making its conjugate a constant of the 
motion. Since the transformation is a near-identity one, the lowest order of 
H* has the same form as that of H , giving 

H*© = z, ♦ £ d k H*«(j[) 

k*J. 


( 42 ) 
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(k) *(k) 

Methods then exist for deriving (J and H ' ' order by order 
Suppose now that two near-identity canonical transformations are given 

* * y + b ' 1 * 

3=1 

» = y + 

8=1 

either of which eliminates the angle variable from the Hamiltonian. Let 
furthermore 

£^(y) * 4^ B \y) (44) 

for 

s = 1 ( 2, ... (k-l) 

Then it will be show^that differs from at most by a conjugate 

gradient of a function (y) independent of y n . 

Clearly such a property would allow the two methods discussed in this 

work to be matched order by order. If w , for instance, represents a solution 

of a given problem by the Poincare-Von Zeipel method and z represents a solution 

by the Kruskal method, then z can be made equal to w to any desired order. 

(k) 

To achieve this one must first choose so as to make jz canonical and 

then add to each order the appropriate £ (y) which makes the 
(k) (k) 

corresponding ^ and ^ equal to each other. Actually, one could 
also work in the opposite direction, since the Poincare-Von Zeipel method also 
con ta i ns a certain amount of arbitrariness in each order, but we shall not 
consider this possibility here. 
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Let 


and 



8=0 




8=0 



(£) 


be the two alternative forma of the Hamiltonian. Sinoe the transformation 
is time-independent 

H *(£) = H**(v) » H(y) (45) 


Also, since it is a near-identity transformation, the 0(l) parte of the above 
equations must be equal, which leads to 




= * 1 


Ve shall furthermore assume (and later justify) 



(46) 


H *(s) = R **(s) 


(47) 


(s = 1, 2, ... k-1 ) 

Substituting the expansions in (45) 

(m) (m) / » 

By means of expansion operators T and 8 (see Z ) this may be 
broken up into a series of equations, one for each order in 51 • The equation 
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for 0(^ ) then is 
k 


> T (m) . H* (kH,) (jf) 


<1 

m=0 


£ S M . H** (k_n) (y) (49) 


B=0 


where * narks "operates on" and where 


r«» = a (0)- 


.( 1 ) 


f (l) - V 


s< l > . 


and in general 


T (b) = + S (a) (|) 

S (b) = v/ b) V - s M (p 


,(») 


(50) 


J 


with N an operator depending only on orders of its argument that are lower 
than the m-th. Substitution and use of (AS) give 

H* (k) + 5 (k > vyi + jif <k) (J). y x - ZL t (, U* M 


m=l 


h"W* ♦ *ZKh“W| 


m=l 


J 

(51) 


The two expressions in curly brackets depend only on lower orders; they are 
therefore equal to each other and nay be dropped, leaving 


h‘« (? ) - H**< k >(5) . (^< k > - ^-vy. 


(52) 
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Now by (24), scalars 



and 



must exist such that 


s" ■ 

. + 

vx (k) 

* (k) - 


vt (k) 


(53) 


tWO ( k ) 

Since lower orders of the^ expansions are equal, the two f vectors are equal 
too, leaving 


£ (k) . y (k) 


y(%^- =■ (54) 


To prove our assertion we must show that ^ is independent of y n , Substituting 
in (52) 


H* (k) (y) - H** (k) (?) = ^4/5^ = Di/?y n (55) 

Now ^ is allows 1 to depend on y n only in a periodic manner, from which it 
follows that ^)$A>y n also depends periodically on y n . However, the 
left-hand side is independent of that variable , so that must be independent 
of y n . This proves the main assertion. Incidentally, one also finds. 

H*< k > = H**( k > (56) 


which justifies equation (47) 
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APPENDIX : THE INVERSE TRANSFORMATION 


Adding ( 5 ) and (24) and cancelling zeroth order terms gives 

2 ^ 'i (k) u) - -2y$ tk) w 

= - y £ k < + z^* \ w <*> ) u-i) 

If V is in z-space and * deflates operation 

$ 00 (**Z t "'l <,,){ *i) = $ <k) w 

= 21 £ a sK 2 k) (i) 


(A-2) 


where 


s‘°> = 1 

.( 2 ) _ 


m=0 




(a-?) 


+ i W 


and so forth. If this is substituted in (A-l), orders of £ may be individually 
equated since everything is now expressed in z . Separating the m = 0 term 
from the rest then brings the 0( £ ) relation to the form 


n (k> ■ 

- ! (k) 

- L s u) * j (kH,) 


m=l 

If all lower orders of 

s w — 

(in) 

known, those of **[ 


(A-4) 


used for constructing the S 


(m) 
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